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We describe structure of the T-matrices, scattering matrices, and Green functions on un- 
physical energy sheets in multichannel scattering problems with binary channels and in the 
three-body problem. Based on the explicit representations obtained for the values of T- 
and S'-matrices on the unphysical sheets, we prove that the resonances belonging to an un- 
physical sheet are just those energies where the correspondingly truncated scattering matrix, 
taken in the physical sheet, has eigenvalue zero. We show, in addition, that eigenvectors of 
the truncated scattering matrix associated with its zero eigenvalue are formed of the breakup 
amplitudes for the respective resonant states. 



I. INTRODUCTION 

Resonances of multichannel systems play a crucial role in various problems of nuclear, atomic, 
and molecular physics. In a wider sense, resonances represent one of the most interesting and in- 
triguing phenomena observed in scattering processes, and not only in quantum physics but also in 
optics, acoustics, radiophysics, mechanics of continua etc. Literature on resonances is enormous 
and in this short introduction we have a chance to mention only several key points in the his- 
tory of the subject and to refer only to a few key approaches to quantum-mechanical resonances, 
necessarily leaving many others a part. 

With a resonance of a quantum system one usually associates an unstable state that only ex- 
ists during a certain time. The original idea of interpreting resonances in quantum mechanics as 
complex poles of the scattering amplitude (and hence, as those of the scattering matrix) goes back 
to G. Gamov [8]. For radially symmetric potentials, the interpretation of two-body resonances as 
poles of the analytic continuation of the scattering matrix has been entirely elaborated in terms of 
the Jost functions [|l2ll . Beginning with E. C. Titchmarsh ll34ll it was also realized that the S-matrix 
resonances may show up as poles of the analytically continued Green functions. 

Another, somewhat distinct approach to resonances is known as the complex scaling (or com- 
plex rotation) method. The complex scaling makes it possible to rotate the continuous spectrum 
of the A^-body Hamiltonian in such a way that resonances in certain sectors of the complex energy 
plane turn into usual eigenvalues of the scaled Hamiltonian. In physics literature the origins of 
such an approach are traced back at least to C. Lovelace [190. A rigorous approval of the complex 
scaling method has been done by E. Balslev and J. M. Combes A link between the 5-matrix in- 
terpretation of resonances and its complex rotation counterpart was established by G. A. Hagedorn 
|@] who has proven that for a reasonable class of quickly decreasing potentials at least a part of the 
scaling resonances for an A'^-body system {N < 4) turns to be also the scattering matrix resonances. 
We remark that the complex scaling seems to be the most popular approach to practical calculation 
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of resonances, particularly in atomic and molecular systems (see, e.g., Refs. IS LZL lilL lla, I2JJ] and 
references cited therein). 

If support of the interaction is compact, the resonances of a two-body system can be treated 



within the approach created by P. Lax and R.Phillips 11811 . An advantage of the Lax-Phillips 
approach is in the opportunity of giving an elegant operator interpretation of resonances. The two- 
body resonances show up as the discrete spectrum of a dissipative operator which is the generator 
of the compressed evolution semigroup. An operator interpretation of resonances in multichannel 
systems, based on a 2 x 2 operator matrix representation of a rather generic Hamiltonian, can be 
found in [20Q. 

For more details on the history of the subject and other approaches to resonances, as well 
as for the bibliography we refer to books [|Il,j3|, 0, |28l 31, 32] (it might also be useful to look 



through the review parts of papers fl^ and psi]). Here we only notice that, in contrast to the 
"normal" bound and scattering states, the resonant ones still remain a quite mysterious object and 
many questions related to resonances are still unanswered. This is partly related to the fact that, 
unlike the "normal" spectrum, resonances are not a unitary invariant of a self-adjoint (Hermitian) 
operator. Moreover, following to J. S. Rowland fiol] and B. Simon f3^, one should conclude that 
no satisfactory definition of resonance can rely on a single operator on an abstract Hilbert space 
and always an extra structure is necessary. Say, an unperturbed dynamics (in quantum scattering 
theory) or geometric setup (in acoustical or optical problems). Resonances are as relative as the 
scattering matrix is itself. 

In our approach present in this report we follow the typical setup where the resonances arising 
due to an interaction V are considered relative to the unperturbed dynamics described by the kinetic 
energy operator Hq. The resolvent G{z) = {H — z)^^ of the total Hamiltonian H = Hq + V is an 
analytic operator- valued function of z G C \a{H). The spectrum g{H) of // is a natural boundary 
for holomorphy domain of G{z) considered as an operator- valued function. However the kernel 
G(-,-,z) may admit analytic continuation through the continuous spectrum of H. Or the form 
{G{z)(p, V^) may do this for any 9, i// of a dense subset of the Hilbert space Sj. Or the "augmented" 
resolvent PG{z)P admits such a continuation for P the orthogonal projection onto a subspace of 
S). In any of these cases one deals with the Riemann surface of an analytical function. 

In the simplest example with H = Hq = —A, the two-body kinetic energy operator in coordinate 
representation, we have 

G{x,x',z) 



Alt \x-x'\ 

where x,x' are three-dimensional vectors. Clearly, G{x,x' ,z) as a function of the energy z has a 
two-sheeted Riemann surface which simply coincides with that of the function z^/^. 

In this way one arrives at the concept of the unphysical energy sheet(s). The copy of the com- 
plex energy plane where the resolvent G{z) is considered initially as an operator- valued function 
is called the physical sheet. The remainder of the Riemann surface is assumed to consist of the 
unphysical sheets (in general, an unphysical sheet may only be a small part of the complex plane). 

Meanwhile, any analytic function is uniquely defined by its values given for an infinite set of 
points belonging to its initial domain and having at list one limiting point. Usually one knows the 
T-matrix or Green function on the whole physical sheet which means that, at least in principle, it 
should be possible to express their values on unphysical sheets through the ones on the physical 
sheet. 

In sm (see also [ l24l I27I1 ) we have found just such expressions. More precisely, we have 



derived explicit representations for the values of the two- and three-body G(z), T{z), and S{z) on 
unphysical energy sheets in terms of these quantities themselves only taken on the physical sheet. 
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The same has been also done for analogous objects in multichannel scattering problems with 
binary channels [id]. The representations obtained not only disclose the structure T-matrices, 
scattering matrices, and Green functions on unphysical energy sheets but they also show which 
blocks of the scattering matrix taken in the physical sheet are "responsible" for resonances on 
a certain unphysical sheet. This result paves the way to developing new methods for practical 
calculation of resonances in concrete multichannel systems and, in particular, in the three-body 
ones (see, e.g. [IjilflEs']). As a matter of fact we reduce all the study of resonances to a work 
completely on the physical sheet. 

The present report essentially extends the presentation given recently in ||23|]. 



II. TWO-BODY PROBLEM 

In general, we assume that the interaction potential v falls off in coordinate space not slower 
than exponentially. When studying resonances of a two-body system with such an interaction 
one can employ equally well both coordinate and momentum representations. However in the 
three-body case it is much easier for us to work in the momentum space (for an explanation see 
iisll . p. 149). This is one of the reasons why we proceed in the same way in the two-body case. 
Thus, for the two-body kinetic energy operator ho we set {hof){k) = k^f{k) where ^ G stands 
for the reduced relative momentum. In case of a local potential we have v{k, k') = v{k — k') and 
v{k)=^i^. 

Surely, we need to add some requirements on the analyticity of the potentials in their complex 
momentum variables as well as on their fall-off as the real parts of the momenta approach infinity 
(see [25] and [26] for details). For simplicity, through all this presentation we assume that the po- 
tentials involved are holomorphic functions of the momenta on the corresponding whole complex 
spaces (that is, on the whole C-^ in the two-body case). 

The transition operator (T-matrix) reads 

t{z)=v-vgiz)v, (2.1) 

where g{z) = (h — z)^^ denotes the resolvent of the perturbed Hamiltonian h = ho + v. The operator 
t is the solution of the Lippmann-Schwinger equation 

t{z)=v-vgo{z)tiz), (2.2) 

that is, in terms of the its kernel we have 

t{k,k,z)=v{k,k)- dq ^ 

q^-z 

taking into account that the free Green function reads 

5{k~k') 



go{k,k',z) 



Clearly, all dependence of ? on z is determined by the integral term on the right-hand side of 
(12.31) that looks like a particular case of the Cauchy type integral 

^{z) = [ dq ■ /^f (2.4) 
Jrn A+q^-z 
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for N = 3. Cauchy type integrals of the same form but for various we will also have below when 
considering a multichannel problem with binary channels in Sec. |in]and the three-body problem 
in Sec. lYl 

Let yix, X G C, be the Riemann surface of the function 

[log(z-A) ifA^iseven. 

If A'^ is odd, is formed of two sheets of the complex plane. One of them, where (z — A)^/^ 

coincides with the arithmetic square root \Jz—'k, we denote by rio. The other one, where [z — 

A) = -y/^- A, is denoted by Di. 

If N is even, the number of sheets of is infinite. In this case as the index £ of a sheet Tlf we 

take the branch number of the function log(z — A) picked up from the representation log(z — A) = 

log|z- A| +i2;r£ + i^ with ^ G [0,2;r). 

Usually the point A is called the branching point of the Riemann surface IH^ • 

The following statement can be easily proven by applying the residue theorem (if necessary, 

consult [|26ll for a proof). 

Lemma 1. For a holomorphic f{g), q G C^, the function 4>(z) given by (|2.4I) is holomorphic on 
C \ [A, +oo) and admits the analytic continuation onto the unphysical sheets Hi of the Riemann 
surface IH^ as follows 

^iz\n,)=^iz)-ini{v^f-^ f dqfi^/^q), (2.6) 

where S^^ ^ denotes the unit sphere in centered at the origin. 

Notice that in (12.61 ) and further on the writing means that position of z is taken on the un- 
physical sheet Uf . If the reference to Ug is not present and we write simply z than one understands 
that we deal with exactly the same energy point but lying on (dropped onto) the physical sheet Ilo. 

Now return to the two-body problem and set 



{go{z)fij2) = dq 
where /i and /2 are holomorphic. Then by Lemma 1 



Mq)f2{q) 
q^-z 



{go{z\ni)fij2) = {gQ{z\uo)fij2)-Tti^z j^^dqfi{^q)f2{^q), 

which means that the continuation of the free Green function go(z) onto the unphysical sheet 111 
can be written in short form as 

gQ{z\u,) =gQ{z)+^Q{z)i\z)i{z), (2.7) 

where ao(z) = —niy/z and j{z) is the operator forcing a (holomorphic) function / to set onto the 

energy shell, i.e. {j{z)f) (k) = f{^k). 

Taking into account (12.71) . on the unphysical sheet Hi the Lippmann-Schwinger equation (12.21) 
turns into 

t' = v-v{go + aofj)t\ t' = t\u,. 
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Hence (/ + v^o)^' = v — Siofjt'. Invert / + v^o by using the fact that t{z) = v — vgot and, hence, 
{I + vgo)-h = t: 

t' = t~aotfjt'. (2.8) 
Apply j{z) to both sides of (|2.8b and obtain jt' = jt — slq jt f jt' , which means 

iT+aojtj^)jt' = jt, (2.9) 

where / stands for the identity operator in ^2(5^) . Then observe that / + ao jt p is nothing but the 
two-body scattering matrix s{z) since the kernel of the latter for z G Ho is known to read 

s{k,k ,z) = d{k,k) - 7ti^/zt{^/zk, y/zk 

Hence 

jt' = [s{z)]-'jt. (2.10) 
Now go back to (|Z8] ) and by using (12.101 ) get t' = t-aot p[s{z)]^^jt, that is, 

t{z\n,)=t{z)-aoiz)tiz)jHzMz)]-'jiz)tiz). (2.11) 

All entries on the right-hand side of (12.1 II) are on the physical sheet. This is just the representation 
for the two-body T-matrix on the unphysical sheet we looked for. 

From (12.1 II ) one immediately derives representations for the continued resolvent, 

^(z|n,) =g + aoiI-8v)p[siz)]-'jiI-vg), (2.12) 

and continued scattering matrix, 

^{z\n,) = ^[s{z)]-'c^, (2.13) 

where S' is the inversion, {S'f){k) = f{—k). Hence, the resonances are nothing but zeros of the 
scattering matrix s{z) in the physical sheet. That is, the energy z on the unphysical sheet Hi is a 
resonance if and only if there is a non-zero vector jz/ of L2{S^) such that 

s{z)£^ = (2.14) 

for the same z on the physical sheet. 

We remark that this fact is rather well known for the partial-wave Schrodinger equations in case 
of centrally-symmetric potentials. For this case the statement that the resonances correspond to 
zeros of the partial-wave scattering matrix si on the physical sheet of the complex energy plane 
follows from its representation (see, e.g., [IJ) 



Slip) = {-I) 



I flip) 

M-p) 



in terms of the lost function fi{p) where / stands for the angular momentum and p for the (scalar) 
complex momentum. This property of si{p) was explicitly noticed in the review article 1122 . 
p. 1357]. Generalizations of the statement to the case of multichannel problems with binary chan- 
nels and to the three-body problem have been given in [26] and [|25l] . respectively. We will discuss 
them below in Sec. ImjandHVl 
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The eigenfunction =2/ in (|2.14l) represents the breakup amplitude of an unstable state associated 
with the resonance z- This means that in coordinate space the corresponding "Gamov vector", i.e. 
the resonance solution to the Schrodinger equation, has the following asymptotics 



exp(iz^/2 


Hi 


^1) 









(2.15) 



-iyfz\x\ 



x = x \x\. 



This claim is a particular case of the statement of Lemma 2 below. 

It should be stressed that the asymptotics (12.151 ) contains no term with the incoming spherical 
wave 

exp(-izi/2|^Jj:|) 

We conclude the section with a remark that in (see also fid. Section 2] and ['22I, Secti- 
on 3]) Yu. V. Orlov was very close to obtaining a representation that would be a version of the 
representation (12.1 II) for partial-wave two-body T-matrices in the case of centrally- symmetric 
potentials. As a matter of fact, only the last step has not been done in [29, 30], the one analogous 
to the transition from equation (12.81 ) to equation (12.111) by using relation (12.101 ). 



III. MULTICHANNEL PROBLEM WITH BINARY CHANNELS 

From now on assume that his anmxm matrix Schrodinger operator of the form 



h 



V21 



V12 



Vim 
Vim 

,('") 



written in the momentum representation. Thus, we assume that 
(a) 4- \n. \ 1.2 



(3.1) 



1,2,, 



,m. 



We restrict ourselves to the case where the channel dimensions satisfy inequahties «a > 3, a = 
1 , 2, . . . , m, and 1 < m < 0°. For simplicity we assume that the potential/coupling terms v^p {ka, k'p ) 
are holomorphic functions of their variables ka G C"" and k'p E , sufficiently rapidly decreasing 
as Reka ^ 0° or Rek'p 0° (see [id]). The thresholds Ai, A2, ■ • • , G M are supposed to be 



distinct and arranged in ascending order: Ai < A2 < 
We also introduce the notations 



< A,„. 





V 







(2) 








\ 



and V : 



/ Vii Vi2 

V21 V22 

\ Vml Vm2 



Vim 
V2m 
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for the unperturbed Hamikonian and the total interaction, respectively. By go{z) and g{z) we 
denote the corresponding resolvents, 



go{z) = {ho~z)-' and g{z) = {h-z'-' 

Like in the one-channel (i.e. two-body) case of Sec. HIl we again begin with the study of the 
T-matrix 

t{z) =v-vg{z)v 
that is the solution to the Lippman-Schwinger equation 

t{z)=v~vgo{z)t{z). (3.2) 

Kernels tap{ka,k'p,z) of the block entries ?c</3(z) of the operator matrix t{z) solve the equation 
system 

^im„y Ay + q^-z 

Similarly to the two-body T-matrix in equation (12.31) . all dependence of the kernels t^p {k, k' ,z) on 
z is determined by the integrals on the right-hand-side of (13.31) . which are the Cauchy type integrals 
just of the form (l24l) . 

In contrast to the two-body case, for m > 2 we arrive at a multi-sheeted Riemann surface with 
number of sheets greater than two. The reason is simple: every threshold Xa, oc = l,2,...,m, 
turns into a branching point. If all the channel dimensions na are odd, the number of sheets should 
be equal to 2"\ that is, in addition to the physical sheet the Riemann surface will contain 2™ — 1 
unphysical ones. If at least one of n^s is even, we will have a logarithmic branching point and 
the number of unphysical sheets will be necessarily infinite. In fact, this Riemann surface simply 
coincides with the Riemann surface of the vector- valued function 

C(z) = (Ci(^),C2(z),...,Uz)), 



where (cf. formula (12.51) ) 



(z-XaY^^ if Ha is odd, 
log(z-Aa) if Ha is even, 



a = 1, 2, . . . ,m. 

To enumerate the sheets of it is natural to use a multi-index 



where each coincides with the branch number for the corresponding function a = 1,2, ...,m. 
In particular, if na is odd then ia may get only two values: either or 1. For even the value of 
ia is allowed to be any integer. The sheets of are denoted by 11^. The physical sheet corresponds 
to the case where all components of i are equal to zero and thus it is denoted simply by Do. 

Each sheet Hi is a copy of the complex plane C' cut along the ray [Ai,+oo). The sheets are 
pasted to each other in a suitable way along edges of the cut segments between neighboring points 
in the set of the thresholds Xa, OC = 1 , 2, . . . , m. In particular, if coming from the sheet n^^^ ^^, 
the energy z crosses the interval A^+i), oc = 1,2, ... ,m, Am+i = in the upward direc- 
tion (i.e. passes from the region Imz < to the region Imz > 0), then it arrives at the sheet 
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n 



I with all indices beginning from ia+i remaining the same while the first a 
indices ij, I < j < oc, change by unity. If nj is odd then i'- = 1 for ij = and £'j = for ij = I; 
if Hj is even then i'j = ij + l. In the case where the energy z passes the same interval {^a^^a+i) 
downward, it arrives at the sheet U/pi oi pi p , p \ where for odd n ; the indices are the same 
as in the previous case and for even rij they change according to the rule 1'^ = ij — I. The indices 
ij with numbers j > oc + l remain unchanged. 

Under the assumption that the kernels tg,p{\/z — Xak,k\z) admit the analytic continuation in 
z through the cuts (the existence of such a continuation may be rigorously approved, see ll26ll ') 
one can perform analytic continuation of the Lippman-Schwinger equation (13.31) from the physical 
sheet Tlo onto any unphysical sheet Yl( of the surface 9^. Of course, the trajectory along which we 
pull z should avoid the branching points Xa- Applying after each crossing the interval (Ai, +oo) 
the corresponding variant of formula (12.61 ) we arrive at the following result 



tap{k.k,z\^)=Vap{k,k)-l^ dq 

y=l>/K„y Ay + q —z 



(3.4) 



I 

7=1 



(^) /„y_i Vay{ka, ^Jz- Xyq) typ ( ^z-Xyq, k'p , z| 



where 



-nn 



Xyfy- 



(3.5) 



Notice that the second integral term on the right-hand side of (13.41) includes the half-on-shell values 
tyfi{\/z — Xyq, ^/3,z|j|P of the T-matrix kernels ?y|3(^,^'^,z|j^ ) taken on the unphysical sheet 11^. 
Thus, like in the two-body case of Sec. |Ill it is convenient to introduce operators jy(z) forcing a 
holomorphic function f{q), q E C'\ to set onto the corresponding energy shell, i.e. 



{jYiz)f){q)=f{^z-Xyq), 7=1,2, 
From these operators we construct a block diagonal matrix 



, m. 



[ h{z) 
h{z) 



\ 







\ 


•y^{z) J 



Using this notation one easily rewrites equation (13.41) in the matrix form 

^Anp)=^-^SQ{z)t{z\^) -vj\z)LA{z)J{z)t{z\^), 
where L and A(z) are diagonal mx m matrices with scalar entries. 



(3.6) 



f ii 

i2 



\ 




\ 



-On 



and A{z) 



fMz) 

A2(z) 



J 



\ 




V 







(3.7) 



A,niz) J 



and J^{z) is the "transpose" of J{z) which means that the product t{z)J\z) has half-on-shell ker- 
nels of the form v^^g {k, Jz — X^k ) . 
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When rearranging (13.61) we first transfer the term vgo(z)K^|nP '-^^ left-hand side of (13.61) and 
obtain 

{l + vgo{z))t{z\^J=v-vjHz)LA{z)J{z)t{z\^^). (3.8) 

Our next step is to invert the operator {l + vgQ{z)) (of course, this is only possible for z not 
belonging to the discrete spectrum of h). Here, we keep in mind that the energy z in this operator 
is from the physical sheet where the Lippmann-Schwinger equation (13.21) holds and thus (/ + 

vgoiz)) = t{z). Using this inversion formula we then derive from (13.81) that 

~tiz)jHz)LA{z)Jiz)t{z\^^). (3.9) 
At this point it is convenient to introduce another diagonal scalar mx m matrix 



L 



/ ii ... \ 

i2 ... 

V ...I,nJ 



(3.10) 



whose diagonal entries are 




if£a = 0, 
if^«^0. 



Clearly, the matrices L, L, and A(z) commute. Moreover, LL = L. Using these facts one rewrites 



(13.91 ) in a slightly different form 

^(4,) -t{z)jHz)LA{z)LJ{z)t{z\^^). (3.11) 

which means that the value Kz|n ) °f '^he T-matrix ? at a point z on the unphysical sheet is 
expressed through the value of t itself taken at the same point z on the physical sheet as well as 
through the half-on-shell value 7(z)?(z|pj ) taken still for z\yi^, and, in addition, multiplied by L 

from the left. Applying the product LJ{z) to both side of (|3.11l) we arrive at a closed equation for 

LJiz)t{z\j,J: 

^+LJ{z)t{z)j\z)LA{z)\LJ{z)t{z\^)=LJ{z)t{z), (3.12) 
where / denotes the identity operator in the sum Hilbert space 

e = L2(5"i-i) ©L2(5"2-i) _ _ . ©L2(5"'"-i). (3.13) 

Therefore, at any point z in the physical sheet where the operator 

Si{z)=f+LJ{z)t{z)j\z)LA{z) (3.14) 

is invertible, we will have 

LJ{z)t{z\^) = [si{z)]-'U{z)t{z). (3.15) 
Notice that si{z) commutes with L, i.e. 

Lsi{z)=Si{z)L, 
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and hence 

LA{z)se{z)-'L = LA{z)se{z)-' . (3.16) 
Taking into account equalities (|3.15l) and (13.161) we obtain from (13.121) the following result: 

t{z\u,) =^iz) -t{z)jHz)LA{z)[se{z)r'LJ{z)t{z) (3.17) 
^t{z)~t{z)j\z)lA{z)Hz)]-\z)J{z)t{z). (3.18) 

These are just the representations for ) we look for: in (13.171 ) and (13.181) values of the 
multichannel T-matrix on an arbitrarily chosen unphysical energy sheet are explicitly written 
in terms of the entries whose values are taken from the physical sheet. Formulas (13.171) and (|3.18l) 
are just the ones that represent a generalization of the two-body representation (12.111) to the case 
of multichannel Schrodinger operators with binary channels. A slightly different version of the 
representations (13.171) and (13.181) was first published in II26I1 . 

The operator matrix S£{z) given by (13.141) is closely related to the total scattering matrix for the 
problem which reads 

s{z)=T+J{z)t{z)j\z)A{z), (3.19) 

Of course, the total scattering matrix contains neither entry L nor entry L. For the matrix S£{z) these 
entries play an important role. Depending on the unphysical sheet Hi under consideration, certain 
rows and columns of the difference matrix [s£{z) — /) = LJ(z)t{z)J\z)L completely consist of 
zero entries. Nullification takes place for those rows and columns of the difference matrix [s{z) — 

l) — J{z)t{z)J\z) whose numbers a are such that the corresponding indices equal zero. This 
is a reason why we call S£{z) the truncated scattering matrix associated with the unphysical sheet 
^^ 

Notice that if instead of (13.21) we start with the transposed Lippmann-Schwinger equation 

t{z) =v-t{z)go{z)v, 

then in the same way we obtain for ) another representation that can be considered as a 
transposed version of the representation (13.171) : 

= -tiz)jHz)L[sliz)]-'Aiz)Uiz)tiz) (3.20) 
^t{z)-t{z)JHz)[sl{z)]-'A{z)U{z)t{z), (3.21) 

where 

sl{z)=T+LA{z)J{z)t{z)jHz)L. 
The operator s^iz) represents the result of truncation of the transposed 5-matrix 

s\z)=T+A{z)J{z)t{z)jHz)- 
From the uniqueness of the analytic continuation by (13.171) and (13.201) it immediately follows that 

t{z)JHz)LA{z)s,{z)-'j{z)t{z) = t{z)JHz)[sl{z)]-'A{z)LJ{z)t{z). 

To describe structure of the scattering matrices s{z) or s^{z) analytically continued to an un- 
physical sheet Hg we need some more notations. First, introduce a block diagonal operator matrix 
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£'{£) of the form S" = diag (^i, ^2, ■ • • , ^m) where S'a is the identity operator on L2(S"«~^) if is 

even and is the inversion, {(^af) (^) = f{~^)^ if is odd. Second, let e(£) be a scalar diagonal 
matrix, e = diag (ei , e2, . . . , e„j) , with the main diagonal entries e^ defined by 

+ 1 forany = 0, ±1,±2,... if is even, 
+ 1 if is odd and la = 0, 
— 1 if Ha is odd and ia = I- 

That is, &a only depend on the corresponding na and ia- It is obvious that if a matrix- valued 
function A (z) is defined on the physical sheet of the Riemann surface 91 by formulas (|3.5I) and 
(13.71) . then after the analytic continuation to the sheet II^ it acquires the form 



(3.22) 



Now we are ready to present our main result concerning the 5-matrices. We claim that after 
continuation to the sheet 11^ their values are expressed by the formulas 



^(^In,) 



1+ t Ae- t LAsn ^tAQ 



/ + qAT - qAT [s\] ^AL V 



(3.23) 
(3.24) 



where we use another shorthand notation 

?{z)=jmz)j\z). 

The argument z of the operator- valued functions S(^{z), s\{z), J{z), ^^(z), and A{z) on the right- 
hand sides of (13.231) (13.241) is a point on the physical sheet Do having just the same position on 
the complex plane as the point zL on the sheet 11^ on the left-hand sides of (13.231) and (13.241 ), 
respectively. 

At last, we present the representation for the continued resolvent on the sheet Hf 



S{z\n,) =S+ {l-gv)PALsj'j{l-vg), 

=g+ {l~gv)j^sl]-'ALj{l-vg). 



(3.25) 
(3.26) 



In this report we skip derivation of the representations (I3.23l) - (|3-26l) . The interested reader 
may find it in 11271 Sections 1 .4 and 1 .5] (see also n26i\ ). Here we only remark that the derivation is 
rather straightforward being based directly on the representations (13.201 ) or (13.211 ) for the T-matrix. 

The most important consequence of the representations (|3.23I) - (I3.26I ) is the fact that all energy 
singularities of the T-matrix, scattering matrices, and resolvent on an unphysical sheet 11^, dif- 
fering of those in the physical sheet, are just the singularities of the inverse truncated scattering 
matrix [i'f (z)]^^ (or, and this is the same, the ones of its transpose [^^(z)]^^). This means that 

resonances on sheet Tli correspond exactly to the points z 
on the physical sheet where the operator S£{z) has eigenvalue zero, 



(R) 



i.e. the resonances on 11/ are those energies z on ITq where equation 



Si{z)£^ = 



(3.27) 
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has a non-trivial solution £/ ^O'm the sum Hilbert space (5 given by (|3.13l) . 

Eigenvectors of the truncated scattering matrices S£{z) associated with resonances have a quite 
transparent physical meaning. Assume that z is a resonance on the unphysical sheet Uf. This 
implies that for the same energy z on the physical sheet Uq equation (13.271) has a solution £/ j^O, 
£^ = {si\,s^2, . . .,s^rn)^. Clearly, the components s^a of the vector .si are non-zero only for the 
channels a such that 7^ 0. Taking into account that (13.271) can be written in the equivalent form 

£^ = -LJt{z)J^LA{z)£/, (3.28) 

this conclusion follows from equality 

(/-L)^ = 0. 

Notice that the latter holds since L{I — L) =0. 

Along with the vector we also consider an "extended" vector that is obtained of ,2/ as a 
result of replacing the projection L on the right-hand side of (13.281) with the identity operator, i.e. 

^ = -Jt {z)J^LA {z)£^. (3 .29) 

Clearly, = Li?r _ _ _ 

We claim that up to scalar factors the components £/\{k\), £^2{ ^2), • • Mn{km) of the eigen- 
vector ^ make sense of the breakup amplitudes of the corresponding resonance state in channels 
1, 2,..., and m, respectively. In particular, these amplitudes determine angular dependence of 
coefficients at the spherical waves in the asymptotics of the channel components of the resonant 
solution to the Schrodinger equation in coordinate representation. 

To give some details, let us denote by and v* the coordinate-space version (Fourier transform) 
of the operators ho and v, respectively. Namely, let 

hQ = diag( Ai - A;cj , A2 - A;^^ ' • • • ' - ^x„, ) , 

where Aj^^ , a = 1 , 2 . . . , m, stands for the Laplacian in variable Xa E R"" . 

In the statement below we restrict ourselves to the case where absolute values of the unphysical- 
sheet indices corresponding to the even-dimensional channels are less than or equal unity, i.e. we 
assume that if na is even then \Ia\ < 1- Recall that if na is odd then automatically = or = 1- 



Lemma 2. Assume that z is a resonance on an unphysical sheet 11^ with multi-index i ~ 
(£i,£25 • • • Am) s^ch that l^al < 1 for all a = 1,2, ... ,m. Let E & be a non-zero solution to 
equation (|3.27l) for the same energy z but belonging to the physical sheet. Then for this z the 
Schrodinger equation 

{h* + v*)w* = zw* (3.30) 

has a non-zero {resonant) solution l//*^, = (l//*^ j , l//*^, 2, • • • , V^res,n)^ whose components V^res,a(-^a) 
for i-ai^^ possess exponentially increasing asymptotics, 

wL,aM = Ca{z,la){^a{-Xa)+o{\)) ^ ^^\m ^ (^.31) 
while for = their asymptotics is exponentially decreasing, 

¥L,aM =^Ca(z,£a)K(Ja)+o(l))p|^^, (3.32) 
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where £^{ka) stand for the corresponding components of the extended vector (13.291) and 

Ca(z,^a) = W2-e ' ''(z-Aa)^ (3.33) 

For the function (z — A^)"^ on the right-hand side of (13.331) mam branch is chosen. 
Complete proof of this statement may be found in ['27', Section 1.6]. 

The functions V^res.a(-^a) taken altogether form the Gamov vector corresponding to the res- 
onance energy z (see, e.g. [4, 28]). Just asymptotic formulas (13.311 ) and (13.321 ) prove that the 

functions s^a{ka)^ 7^ 0^ and ^/a(^a), ta = 0, make sense of the breakup amplitudes describing 
decay of the resonant state along open and closed channels, respectively. 



IV. THREE-BODY PROBLEM 

In this section we give a sketch of our results on the structure of the T-matrix, scattering ma- 
trices, and Green function on unphysical energy sheets in the three-body problem. For detail 



exposition of this material see Refs. Il25ll or ll27ll 



Let be the three-body kinetic energy operator in the center-of-mass system. Assume that 
there are no three-body forces and thus the total interaction reads V = vi -|- V2 + V3 where v«, 
a = 1,2,3, are the corresponding two-body potentials having just the same properties as in Sec. 

m 

The best way to proceed in the three-body case is to work with the Faddeev components 

= 5a/3Va-v«G(z)v^ (a,/3 = 1,2,3) 

of the T-operator T{z) = V — VG(z)V where G(z) denotes the resolvent of the total Hamiltonian 
H = Hq + V . The components Mf^p satisfy the Faddeev equations 

Mapiz) = dapta{z)-ta{z)Go{z) £ M^p{z) (4.1) 
with Go{z) = {Ho-zy^ and 



where ka.Pa denote the corresponding reduced Jacobi momenta (see \\25\\ for the precise definition 
we use) and P = {ka^Pa) ^ is the total momentum. 

Assume that any of the three two-body subsystems has only one bound state with the corre- 
sponding energy Cce < 0, a = 1,2,3. Assume in addition that all of these three binding energies 
are different. It is easy to see that the thresholds £1, £2, £3, and are associated with particular 
Cauchy type integrals in the integral equations (14.11) . By Lemma 1 the two-body thresholds £« 
appear to be square-root branching points while the three-body threshold is the logarithmic one. 
In order to enumerate the unphysical sheets we introduce the multi-index £ = (£0,^1,^27^3) with 
£0 = ...,—1,0,1,... and = 0, 1 if a = 1 , 2, 3 . Clearly, only encircling the two-body thresh- 
olds one arrives at seven unphysical sheets. The three-body threshold generates infinitely many 
unphysical sheets. (There might also be additional branching points on the unphysical sheets, in 
particular due to two-body resonances.) 
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It turns out that the analytically continued Faddeev equations (14. lb can be explicitly solved in 
terms of the matrix M = {M^^p} itself taken only on the physical sheet, just like in the case of the 
two-body T-matrix in Sec. |ll]and multichannel T-matrix in Sec. [nil The result strongly depends, 
of course, on the unphysical sheet 11^ concerned. More precisely, the resulting representation reads 
as follows ^ 

M\n,=M + QMLSfLQM. (4.2) 

In the particular case we deal with, L and L are 4 x 4 diagonal scalar matrices of the form 
L = diag(£o,^i, ^2,^3) and L = diag(|£o|, ^1,^2,^3)5 respectively; 5^(z) = I + L{S{z) — I)L is a 
truncation of the total scattering matrix S{z) and the entries Qm, Qm are explicitly written in terms 
of the half-on-shell kernels of M (see formula (7.34) of [25]). From (14.21 ) one also derives explicit 
representations for G{z\uf ) and S{z\ut) similar to those of (I3.25|) and (13.231) . respectively. 

Thus, to find resonances on the sheet Tl^ one should simply look for the zeros of the truncated 
scattering matrix Se{z), that is, for the points z on the physical sheet where equation Sf{z)£^ — 
has a non-trivial solution The vector £/ will consist of amplitudes of the resonance state to 
breakup into the various possible channels. Within such an approach one can also find the three- 
body virtual states. 

In order to find the amplitudes involved in S(, one may employ any suitable method, for example 



the one of Refs. [|13l . Il4l 11511 based on the Faddeev differential equations. In these works the 



approach we discuss has been successfully applied to several three-body systems. In particular, 
the mechanism of emerging the Efimov states in the ^He trimer has been studied il3l [l5n. 
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